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1 INTRODUCTION 



ABSTRACT 

We identify migrating stars in an N-body hybrid simulation of a Milky- Way-like disk. 
Outward migration can occur when a star in a low eccentricity orbit lags a short- 
lived local spiral arm density peak. We interpret short lived local density peaks, that 
appear and fade on approximately an orbital period, as arising from positive inter- 
ference between spiral density wave patterns that are longer lived. We find that stars 
near such a peak can migrate over a significant distance in galactocentric radius dur- 
ing the peak lifetime, providing that the peak is sufficiently dense. We propose that 
short lived density peaks, caused by interference between spiral density waves, can 
induce radial migration even when there arc no spiral density waves that are strong 
near their corotation resonance. Using a Gaussian bar model for the potential pertur- 
bation associated with a narrow transient spiral feature, estimates of the migration 
rate, angular offset between particle and spiral feature, and maximum eccentricity for 
migrators roughly agrees with the values measured in our simulation. When multi- 
ple spiral density waves are present, local density peaks can appear and disappear 
on timescales faster than the timescale estimated for growth and decay of individ- 
ual waves and the peak surface density can be larger than for any individual wave. 
Consequently, migration induced by transient density peaks may be more pervasive 
than that mediated by the growth and decay of individual patterns and occurring 
at their corotation resonance. We discuss interpretation of transient-like behavior in 
terms of interfering patterns, including estimating a coherence time for features that 
appear due to constructive interference, their effective angular rotation rates and the 
speed and direction that a density maximum would move across a galaxy inducing a 
localized and traveling burst of star formation. 



Bar and spiral structures can change the mean orbital ra- 
dius or guiding centers of stars and so migrate them from 
their birth radius (|Wielen|1977||Fuchs fc Wielen|1987||Wie- 
len et al.|[T996] [SeTTwood fe Binney||2002| |Roskar et al 



migration, has been modeled as a diffusive stochastic process 
( Jenkins fc Binney||1990] |Schonrich fc Binney1|2009| |BlancT 



2008[ |Minchev et al.||2011| |Roskar et al.||2012| |Grand et 
al. 2012a|b I. The resulting radial mixing is required to ac- 



count for stellar metallicity distributions (ISchonrich & Bin- 



ney 20091. Many stars, including the Sun, (Adams 20101, 



were born in a cluster, (Eggen 1992 De Silva et al. 2007 



Hawthorn et al.|2010 1 . One difficulty of this approach arises 
from the small number of Galactic rotation periods in a Hub- 
ble time. The rotation period of the Sun is about 250 Myr, 
giving only 40 rotation periods during the lifetime of the 
Galactic disk. If spiral arms require a few rotation periods 
to grow and decay, then only about a dozen non-interfering 
patterns can arise during the lifetime of the Galactic disk. 
Recent studies suggest that multiple patterns can co-exist 



Bubar & King 2010 1 that then dissolved or dispersed to 
other radii in the Gal axy fWielen|1977| |Wielen et al.||1996 



in a galactic disk (Henry et al.||2003 


Naoz & Shaviv||2007 


Meidt et al.|2009 Quillen et al.|2011 


Sellwood|2011l. Only 



Portegies Zwartl[2009l |Brown et al.||20ib| |Bland-Hawthorn| 



et al.| 2010 1. Understanding migration and heating mecha- 
nisms in the stellar disk is needed to interpret forthcoming 
surveys of stars in the Galaxy. 



Transient spiral structure (Sellwood & Carlbert 
Toomre||l990l |Toomre fc Kalnajs||1991| |Baba et~aT 



1984 



2009 



if multiple patterns grow and decay in an uncorrelated or 
non-interfering way can stochastic models for radial migra- 
tion be adapted in a straightforward manner to model radial 
migration. Interference between patterns may increase the 
heating rate ( jMinchev fc Quilen|2006| |Minchev et al. 12012] ), 



Fujii et al.||2011| |Sellwood||20U[ ), causing both heating and 



induce chaos in the dynamics ( Quillen||2003 1, account for 
star formation in armlets ( Henry et al.||2003 l, gaps in local 
velocity distributions ( Quillen et al.|2011 1, and modify disk 
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surface brightness profiles ( Minchev et al.||2012 l. We focus 
here on how the presence of multiple spiral density waves 
influences stellar migration. 

Stellar migration can occur when a star is temporarily 
trapped in the corotation region (or resonance) of a spiral 
density wave ( |Lynden Bell fc Kalnajs|1972||Sellwood fc Bin- 
ney||2002 1. We outline the physical mechanism as proposed 
and described by |Sellwood fc Binney| ( |2002[ ). The star is 
captured into the corotation resonance as the spiral density 
wave grows. When there is a fixed rotating spiral pertur- 
bation, the Jacobi integral, Ej — E — LQ. a , is conserved. 
Here E and L are the energy and z component of the star's 
angular momentum per unit mass, and £l 3 is the pattern 
speed of the spiral density wave. The star can vary in angu- 
lar momentum, L, while the spiral density wave is strong. 
After the spiral density wave dissipates, the star can be left 
at a different angular momentum, but may not have sig- 
nificantly increased in eccentricity. A number of numerical 
studies have supported this migration mechanism, includ- 
ing the seminal work |Sellwood fc Binney| ( |2002[ > , but also 



Roskar et al. (2008) using SPH simulations, and recently 
Grand et alf(|2012a|b[l in both N -body and SPH/N-body 
simulations. |Grand et aL (2012a b) illustrated that outward 
migrating stars lagged a spiral density peak and vice-versa 
for inwards migrating stars. They also showed that migrat- 
ing stars remained on nearly circular orbits, as would be 
expected when Ej is conserved and Q m Q a during migra- 
tion. Here Q is the angular rotation rate of the migrating 
star. Distributions of the change in angular momentum as 
a function of initial angular momentum, AL vs Lq, from 



numerical simulations ( 


Sellwood & Binney||2002 Roskar et 


|al.|2012||Minchevet al. 


20121), exhibit features at corotation 



resonances, with positive variations in AL occurring inside 
corotation and vice-versa outside of the resonance. The the- 
ory (as outlined by |Sellwood fc Binney|[2002[ ) for migration 
due to a corotation resonance is simplest in the presence of 
a single slowly growing and then fading spiral pattern. We 
examine here how this mechanism might be modified in the 
presence of multiple spiral density waves. 

In this paper we examine migration in a simulation with 
multiple spiral density waves. We study individual migrating 
particles to probe the mechanism accounting for their rapid 
variations in angular momentum. We then discuss modifica- 
tions to the corotation mechanism for radial migration in a 
setting with multiple spiral density waves. 



2 MIGRATION IN AN N-BODY SIMULATION 

We use the disk galaxy simulation previously described and 



studied by Quillen et al. ( 2011 1. The initial conditions for the 



simulations were for a model Milky Way galaxy, generated 
with numerical phase phase distribution functions using the 
method discussed by Widrow et al. (120081). Spectrograms 



were used to identify a bar and three spiral patterns in this 
simulation, two with two arms, (m = 2), and one with three 
arms (m = 3, see Table 1 by Quillen et al.|20lT I . The spiral 
structure (shown in log r vs angle) shown in Figure 3 by 
IQuillen et al] p0lT| can be compared with similar figures 
of other studies (see Figure 9 by Grand et al. 2012a and 
10 by Grand et al. 2012b I. While m = 2 and 3 structures 




Figure 1. Histograms of the numbers of stars that experienced a 
change in angular momentum, AL (shown on the j/-axis), versus 
initial angular momentum, Lq (x-axis), for the disk particles in 
the simulation. From left to right the histograms are shown every 
0.3 Gyr. Here Lq and AL are in units of 100 km s _1 kpc. 



the one studied by Grand et al. (2012a I is dominated by 



patterns with higher numbers of arms (m between 3-7) and 
has a flocculent rather than grand-design morphology. 

In Figure [l] we show a histogram of changes in angular 
momentum, AL, computed every 0.15 Gyr, as a function of 
initial angular momentum, Lq, in units of 100 km s" 1 kpc for 
the disk particles in the simulation. As have previous stud- 
ies (e.g., Roskar et al.|[2012 Minchev et al.(2012 |, we see a 
strong feature associated with the bar's corotation radius at 
Lq ~ 8 km s _1 kpc, moving progressively outwards as the 
bar slows down. However, changes in angular momentum 
occur all over the galaxy and changes in angular momen- 
tum persist after the bar has grown. This migration cannot 
be mediated by the growth and decay of the bar's corota- 
tion resonance as the bar does not decay. We conclude that 
migration is pervasive and on-going. Because of the num- 
ber of spiral density waves, this is perhaps not surprising. 
However, the spectrograms (see Figures 4-6 by |Quillen et al.| 



2011 1 show that the stronger patterns in the simulation con- 
tain little power at their corotation radii. Only the bar and 
the m — 3 pattern contain power near corotation (see Figure 
6 by Quillen et al. 2011 1. A similar problem was discussed by 



dominate our simulation and that by |Grand et al.| ( |2012b[ ) , 



Grand et al. (2012a I whose spectrograms showed that most 
spectral features were strong at frequencies above or below 
corotation (see their Figure 6). Nevertheless they also found 
that radial migration occurred in their simulation. Grand 
et al.| ( |2012b[ ) illustrated extreme migration near the end of 
their bar, and we see from the AL vs Lq histograms that 
migration in our simulation occurs in the same region. 



2.1 Finding migrators 

We identify particles that migrate large distances in the fol- 
lowing manner. We first identified particles with an over- 
all change in angular momentum of 15 < |AL| < 50 
km s _1 kpc. Of these 'extreme migrators', a randomly se- 
lected sample were inspected for radial movement in prox- 
imity to the bar and spiral arms. Of these we chose a single 
fiducial particle to study in more detail. In Figure [2^l we 
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show the trajectory of a particle that migrated from radius 
r = 4.1 kpc at 500 Myr to 8.6 kpc at 750 Myr, along with 
39 other particles that were near, (within 0.1 kpc in galacto- 
centric radius and 5° in azimuthal angle) the chosen particle 
at 500 Myrs, the start of the chosen particle's outward mi- 
gration. The positions of these particles are shown at 525 
Myr (just after they are all in the same region) and at 25 
Myr intervals afterwards. At 500 Myr, all 40 of these parti- 
cles have similar angular momentum (within 25 km s _1 kpc 
of the chosen particle), however they do not all have the 
same radial velocity. We use u to denote the radial veloc- 
ity component in galactocentric coordinates with convention 
positive u for a particle moving toward the galactic center. 
In Figure^ particles with u ^ —30 km s _1 are shown as 
asterisks and include the chosen particle, the filled squares 
are those with u ^ 30, and the open squares with |u| < 30 
km s -1 . Even though the chosen particle is in a nearly cir- 
cular orbit (and has little epicyclic motion) it has a negative 
u corresponding to an outward drift. 

The positions of the 40 particles are shown in Figure 
[2^i along with the disk density in cylindrical galactocentric 
coordinates. The azimuthally averaged value of the surface 
density has been subtracted at each radius. The rs-axis is the 
azimuthal angle, 9, in degrees and the y-axis is log 10 radius 
in kpc. Rotation is in the positive 9 direction (moving to the 
right in Figure [2^) and the spiral structures are trailing with 
a pitch angle of approximately 24° (as measured previously 
by |Quillen et al.|2011[ ). During the time 500 to 750 Myr the 
chosen particle, and many of the particles with similar initial 
velocities, lag just behind a strong local density feature and 
it is during this time that the particle moves outwards. The 
local density peak exerts a force in the azimuthal direction 
on the particle. This force gives a torque that causes the 
particle's angular momentum to increase. As long as the 
particle remains approximately the same distance away from 
the density peak and lags the peak, it continues to gain 
angular momentum and moves outwards. 

Our fiducial or chosen particle is used in a number of 
Figures. Shown in Figure [3] is the surface density as a func- 
tion of angle, at three different times during the migration 
for our selected particle. The angular difference between this 
particle and nearest density peak is shown in Figure [4] The 
density of the nearest peak and radial migration rate for the 
same particle are plotted in Figure [5] 

At every timestep we measured the angular distance be- 
tween our fiducial particle (one of the asterisks in Figure [2^) 
and the nearest density peak. We denote the angle difference, 
A9 = 8 P eak — dp, between the azimuthal angle of the nearest 
spiral density peak, 9 pea k, and the particle, 9 P ; We measure 
Speak by computing the surface density as a function of az- 
imuthal angle in an annulus with mean radius equivalent to 
that of the particle and with radial width 0.2 kpc. Within 
this annulus we located the highest local maximum disk sur- 
face density that was less than 45° away from the particle, 
and recorded the angle of this peak. Shown in Figure [3] is 
the surface density as a function of 9 — 9 P , at three differ- 
ent times during the migration for our selected particle. At 
each time there is a nearby (within 45°) density peak. The 
particle lags the nearest density peak and this is consistent 
with its outwards migration. Figure [3] also shows the surface 
density in units of Mq pc~ 2 . In subsequent discussion we 
will use the surface density of the nearest peak to estimate 



the migration rate of the particle. The estimated migration 
rate will then be compared to the actual migration rate. 

Plotted in Figure|4^i is the rate of change of angular mo- 
mentum, L, in units of kpc 2 Myr -2 , as a function of time, for 
the selected particle along with the the difference between 
the azimuthal angle of the nearest spiral density peak and 
the particle, A9. On the bottom panel of the plot we also 
show the particle's radius as a function of time. The angular 
rotation in this simulation has positive d9/dt corresponding 
to counter-clockwise rotation. Consequently our plotted an- 
gle, A9, is positive when the particle lags the peak. We see 
from Figure [4^i that the particle moves outwards (L is posi- 
tive and radius increases) when the particle lags the nearest 
density peak (A9 > 0). 

Our illustration that a particle migrates outward when 
it lags a density peak is similar and consistent with the sim- 
ilar illustrations by Grand et al. ( 2012a|b |. The corotation 



resonance model for migration ( Sellwood & Binney 2002 \ 
predicts that the particle would be out of phase with a po- 
tential perturbation while it drifts outwards. In this manner, 
the offset between particle and density peak seen in our sim- 
ulation is consistent with the corotation resonance model for 
migration. 

Figure |4^i also shows the particle radius as a function of 
time. When the particle's eccentricity is non-zero the parti- 
cle oscillates radially (known as the epicycle). Previous stud- 
ies have found that corotation resonances are only effective 
at causing migration among low eccentricity particles (see 
Figure 12 by Minchev et al.|2012 and associated discussion). 
We can see in Figure |4j that the particle only migrates when 
the particle eccentricity is low and that the eccentricity does 
not increase during migration. The positive radial velocity 
value for migrating particles was due to their outward drift 
rather than epicyclic oscillations. The corotation resonance 



model for migration by Sellwood & Binney (2002) is consis- 



tent with the lack of eccentricity variation during migration. 
Figure [2^, shows that particles that have radial velocity u dif- 
fering by 30-40 km s _1 from our extreme migrators fail to 
migrate. For particles nearly in circular orbits, eccentricity 
can be estimated from a maximum u value as e ~ u/v c . 
Hence we crudely estimate that only particles with eccen- 
tricity e < 0.2 migrate in the vicinity of our fiducial particle 
and a maximum migrator eccentricity of e m ax ~ 0.2. 

Figure [4]b we show a particle from Figure [2^, which does 
not migrate significantly. As also seen in the case in which 
the particle migrated, this particle has positive L when lag- 
ging the nearest density peak, but also negative L when 
leading the nearest density peak. As these periods of lag- 
ging and leading occur for similar durations and distances 
to the peak location, this particle does not migrate outward 
like the particle shown in Figure [4^,. 

2.2 The lifetime of a local density peak that can 
cause a star to migrate 

The density feature or peak responsible for the migration 
of the star shown in Figure [2^ appeared at a time of about 
475 Myr and disappeared at a time of about 750 Myr and so 
lasted about 275 Myr. In this simulation the circular velocity 
is about 200 km s _1 and so the rotation period at a Galac- 
tocentric radius of 10 kpc is about 300 Myr. The local peak 
(causing migration of the particle with trajectory shown in 
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525 My r 
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Figure 2. a) Position of 40 stars at different times. These 40 particles were all within 0.1 kpc of each other at 500 Myr, just before the 
top left panel. Each panel shows a different snapshot with snapshots plotted every 25 Myr, starting at 525 Myr. The stars are plotted 
on top of the surface density distribution in cylindrical coordinates. The x-axis gives azimuthal angle and the y-axis logio of the radius 
in kpc. In the plot, the stars marked with an asterisk have a radial velocity, u ^ —30, the filled squares u ^ 30, and the open squares 
\u\ < 30 km s — 1 . The local density peak inducing migration lasted less than 300 Myr. Rotation is with positive dO/dt so patterns move 
to the right. Low eccentricity stars migrate outwards when they lag (and follow) a local density peak in the disk. Only particles with 
small positive radial velocities move outwards. Particles with higher epicycles (the filled and many of the open squares) fail to migrate, 
b) Surface density distribution similar for the same times as a) but for the model described in section [2 . 4| with the parameters listed in 
Table □ 
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Figure 3. Surface density for our selected particle involved in outward migration (one of the asterisks in Figure[2^,), at three separate 
times during migration. The x-axis shows the azimuthal angle, 9 — 8 P , with respect to the angle of the particle. Surface density is 
calculated in an annulus at the particles radius with a width of ±0.1 kpc and is given in units of Mq pc -2 . Rotation is in the positive 9 
direction so the particle lags the nearest density peak for the three times shown. Vertical lines at ±45° show region in which we measure 
the location of the nearest peak and its density. The angular difference A9 between particle and nearest density peak is shown in Figure 
[4] The density of the nearest peak and radial migration rate for the same particle are plotted in Figure [5] 
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Figure |2|i audQJi) existed less than a rotation period. Grand 
et al. ( 2012a I followed a single spiral density peak at a radius 



of about 5 kpc for 120 Myr at which time it disappeared. 
The spiral feature that they traced had a pattern speed (or 
angular rotation rate) that varied with radius and appeared 
to be corotating (see their Figure 5). They demonstrated 
that it induced radial migration in nearby stars. Grand et al. 
( 2012a I commented, "if there are indeed several wave modes 



present, it is evident that they must conspire in a specific 
way in order to produce a spiral arm feature that is ap- 
parently corotating ." Grand et al. (2012a I also commented 



the lifetime of the spiral feature they followed was short. A 
similar short lived spiral feature was traced by |Grand et al.| 
( 2012a I with a lifetime of 160 Myr near the end of a bar. The 



stars identified by Grand et al. (2012b I in a barred galaxy 



simulation with morphology similar to ours also migrated 
when they were near a short lived density feature. 



2.3 



Transient features caused by interference 
between patterns 



Swing amplification spiral models (e.g., Fuchs|20"01 1 require 
a few rotation periods for growth and fading of the spiral 
density wave. However, a density feature that is short lived 
can be produced when two spiral density waves interfere. 
When density peaks are coherent, a large density peak can 
be produced from the sum of the two amplitudes. This den- 
sity peak should survive until the two spiral density waves 
are out of phase. Consider two spiral density waves with m\ 
and m 2 arms each and patterns speeds, Sli and Q2, respec- 
tively. In the frame of the first pattern, the other pattern 
has speed Q2 — fii- If initially both patterns have maximum 
at the same angle, after a time 



t = 



mi 1(^2 — ' 



(1) 



the second pattern would be 180° out of phase at the lo- 
cation of the first pattern peak. We can also consider how 
long it takes for the second pattern peak to reach the next 
minimum of the first pattern. The minimum of these two 
timescales gives a symmetrical coherence timescale 



tr. 



7T 1 

max(mi,m2) |f^2 — fii| 



(2) 



Interference peaks are long lived, (long t co h), only if the pat- 
tern speeds of the interfering waves are similar. For patterns 
with more than one arm, the lifetime of a interference peak 
is shorter than a rotation period at radius ro as long as 



fi 2 - fii 



2max(rai,m2) > 1, 



(3) 



where 0.$ is the angular rotation rate of a particle in a circu- 
lar orbit at ro. We can see if this condition can be satisfied 
in the vicinity of the end of the bar end in our simulation. 
We take pattern speeds measured previously in our simu- 



lation (listed in Table 1 by Quillen et al. 20111. Using a 
bar angular rotation rate Q.b ~ 40 km s~ kpc - and a two 
armed inner pattern with £l 3 « 30 (in the same units) and 
f2o ~ &b(ro/rb) where rb is the bar's corotation radius, we 
estimate 



ilb — 



fit) 



ro 



(4) 



using m — 2 for both patterns. Thus condition equation |3| 
is satisfied outside bar's corotation radius (ro > r;,) for the 
bar and the inner two-armed spiral pattern. The condition 
for a coherence time shorter than a rotation period would 
also be satisfied by a slower three-armed spiral pattern in- 
terfering with the bar. 

The bar in our simulation has a length of 4 -5 kpc 
(corresponding to 0.6 -0.7 in log). However in the simulation 
the bar appears to be longer at time 550 Myr than at 700 
Myr. If the feature that caused migration in our simulation 
is a result of interference between a bar and a spiral density 
wave, then we expect both spiral and bar are in phase during 
the migration. The addition of spiral waves near the end of 
the bar might also account for the some of the variations 
in morphology near the end of the bar, explaining why the 
bar appears to vary in length and how its ends curve at 
different times. At time 500 Myr, the bar extends to 0.8 
(in log). Between 525 and 650 Myr, in the region between 
0.6 and 0.9, one end of the bar progressively curves and 
becomes more tightly wound. Previous studies have also seen 
spiral features that appear to become increasingly tightly 



wound (e.g., Grand et al. 2012b I. In other words the winding 
angle of a density peak appears to tilt as a function of time. 
Another way to describe this behavior is with a measure for 
the angular rotation rate of the density peak that depends 
on radius. A spiral feature that becomes more tightly wound 
also has a density maximum with angular rotation rate that 
varies with radius. Below, we will discuss interpretation of 
spiral features that wind up (increase in winding angle in 
time) or have angular rotation rate that is dependent on 
radius in terms of interference between steady patterns. 

2.4 Describing the simulation density peaks in 
terms of interfering waves 

We can test the possibility that transient structures in our 
simulation are due to spiral density wave interference by see- 
ing if we can mimic the appearance of the density pertur- 
bations, seen in the simulation, with a simple interference 
model. In the same projection as shown in Figure [2^,, we 
have summed three density perturbations each with a fixed 
pattern speed. Each perturbation is described as a surface 
density perturbation with a Gaussian amplitude 

£(r,0,t) = S s exp f- (1 ° Sl °^~ ^ ) /(<MM,*)) (5) 
with function 

f(4>) = cos(0) + 0.25 cos(20). (6) 

For each pattern the parameter X* describes the density 
amplitude. The function, f ((/>), is somewhat more peaked at 
(f> — than a cosine and allows sharp features to be present 
in the surface density. The amplitude peaks at a radius r< 
such that log 10 n = U, consequently U is a parameter that 
describes the mean of the radial region in log space that 
contains the wave. The amplitude extends a radial range 
controlled by the parameter Wi . We define a spiral argument 
for the i-th spiral density wave 



<f>i(r, 9, t) = a, \n(r/r b ) - rrii{0 — Qi(t - t ) - 0,). 



The number of arms or azimuthally symmetry (for a bar) is 
set by rrii. The pattern speed is set by Hi and oti controls the 
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Figure 4. a) For a particle that migrates we plot the rate of 
change of angular momentum (L) as a function of time as a blue 
line with axis on the right hand side. Also plotted, as a red line, is 
the angle to nearest density peak (A9) with axis on the left hand 
side. Note that outwards migration (positive L) only only occurs 
when the particle lags the density (positive A9). We also show 
particle radius as a function of time in the lower panel. Outwards 
migration primarily occurs when radial oscillations are not large. 
Units for L are kpc 2 Myr -2 . b) Similar to a) except for a particle 
that does not migrate a significant amount during the simulation. 



winding angle. The bar perturbation has ai — but when 
a / 0, the spiral perturbations are logarithmic. A logarith- 
mic spiral density wave with rrii arms can approximately be 
described as having surface density and potential dependent 
on the angle <f>i and with amplitude slowly varying with ra- 
dius. The offset of each pattern at time t = to is set by a 
radius, r-,, and angle Bi. Here the winding angle, ji, depends 
on rdO/dr in a density peak with 



cot 7* = ft = 



(8) 



A model that contains two spiral patterns and a bar is 
shown in Figure [2Jd with parameters listed in Table [l] The 
listed parameters are for each pattern: a pattern speed (Hi), 
number of arms (rrii), width parameter (u>i), winding pa- 
rameter (ai), log radius of maximum (k) and angular offset 



(0i). Figure]^), showing the model, can be compared to Fig- 
ure [2^, showing the simulation. We have assumed rotation 
with d6/dt = 0, and trailing spiral patterns corresponding to 
negative ai. The pattern speeds listed in Table [l] are given 
in terms of the bar's pattern speed and the offset angles 
given in terms of the bar's orientation at to = 550 Myr, 
corresponding to the top left panels in Figures [2^, and b. 
The patterns speeds, f2j, are given in units of the bar's an- 
gular rotation rate that is 2.3° Myr -1 (equivalent to the 
0.040 radians Myr -1 or 40 km s -1 kpc -1 measured from a 
spectrogram and listed in Table 1 by Quillen et al. 20111; 
(to convert from radians Myr -1 to km s - kpc - multiply 
by 1000). At early times, (top panels in Figure [2Jd) the spi- 
ral features reasonably well match the simulation with mor- 
phology shown in Figure [2^. However, at later times, the 
model and simulation are increasingly different. The model 
is likely too simplistic. To better match the simulation we 
would probably require a model containing additional pat- 
terns, and with patterns that are not fixed in time but slowly 
vary as the bar slows down. 



The two-armed structure in the model has an angu- 
lar rotation rate of 34.4 and the three armed one has 28 
km s _1 kpc -1 . These patterns were not listed in Table 1 
by |Quillen et al.|2011| as prominent patterns in the simula- 
tion. However they do correspond to fast features seen in the 
spectrograms. The three armed structure corresponds to a 
peak at about 5 kpc at an angular frequency of 0. 090 radians 
Myr - shown in the spectrogram of Figure 6 by Quillen et 



al. (2011 1 (angular frequency is m times the pattern speed). 



The two armed one corresponds to a fainter feature just past 
the end of the bar at an angular frequency of 0.07 Myr -1 
present in the spectrogram shown in see Figure 4 by Quillen 
et al. (2011 1, and that may be associated with a clearer fea- 



ture in the m = 4 spectrogram (their Figure 5) at an angular 
frequency of about 0.14 Myr -1 . 



The reasonable correspondence between the model Fig- 
ure (2jj and the simulation in Figure [2^, the two figures sug- 
gests that an interference model can reproduce transient spi- 
ral features in this simulation. The feature causing migra- 
tion of the star shown in Figures^ (and crudely modeled in 
Figure |2j}) likely appears and disappears due to interference 
between the fast two and three-armed structures with pat- 
terns listed in Table [l] and are consistent with faint peaks 
seen in the spectrograms. 



Using the the bar and two-armed spiral pattern speed 
(and m = 2 for both) listed in Table [I] and equation |2| 
the coherence timescale for peaks arising due to positive in- 
terference is approximately 1.8 bar rotation periods where 
a bar rotation period is approximately 160 Myr. Similarly 
the coherence time between bar and three-armed and two 
and three armed structures is 0.6 and 1 bar rotation peri- 
ods. These short coherence timescales are consistent both 
with the timescale we estimated for the lifetime of the peak 
causing stellar migration and discussed previously in section 
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Table l. Density Wave Model 



Pattern 




mi 




u 


Wi 


8> 


Bar 


1.00 


2 





0.33 


0.23 


0° 


Two- armed 


0.86 


2 


-3 


0.71 


0.14 


-10° 


Three- armed 


0.70 


3 


-6 


0.90 


0.17 


10° 



A three pattern density wave model with parameters chosen to 
match the simulation at around 500 Myr. Here angular rotation 
rates for each pattern, Qi are given in units of the bar's angular 
rotation rate. The number of arms or symmetry in the case of 
the bar is given by m. The winding angle is determined from ai 
and rrii using equation j8j. The density peak is at a radius such 
that li = log 10 r. Each pattern extends over a radial range set 
by Wi . Each pattern is described by Equations [5] - [7] and using 
log 10 r;, = 0.6 with in kpc. The amplitude of each structure 
has = 0.7 in units of the azimuthally averaged value. Here the 
offset angle, 9i, is given in terms of that for the bar that has major 
axis at 130° at time to = 550 Myr, corresponding to the second 
panel in Figure [2k. The resulting density is shown at different 
times in Figure [2p. Parameters were adjusted so that the density 
field is similar to that seen in the simulation and shown in Figure 



EH. 



3 PROPERTIES OF PEAKS ARISING FROM 
INTERFERENCE BETWEEN TWO SPIRAL 
DENSITY WAVES 

The corotation model for migration is consistent with the 
low eccentricities of migrating particles and that migration 
primarily takes place when a particle lags a density peak. 
However our spectrograms don't show patterns (other than 
the bar) that are strong in their corotation region. Further- 
more the peaks that we see in our simulation causing mi- 
gration survive less than rotation period and so are likely 
caused by interference between spiral density waves. Spiral 
patterns need not be long lived, however it is unlikely that 
spiral patterns appear and decay on timescales less than a 
rotation period. We consider the possibility that local peaks, 
caused by positive interference between two patterns, can 
cause stars to migrate. We first consider properties of local 
density peaks produced by interference. We then modify the 
corotation resonance model for migration in this setting. 



3.1 Angular rotation rate of an interference peak 

Using a model for two density waves we now compute the 
angular rotation rate of a local density peak caused by posi- 
tive interference between the two waves. Spiral density waves 
are often described in terms of the number of arms, an am- 
plitude (that is a function of radius), a winding angle and a 
pattern speed. Except for the number of arms, each of these 
parameters could vary in time, and their radial dependence 
could also vary with time. This is a large and poorly con- 
strained multidimensional parameter space for variations. 
To separate between variations caused by single spiral pat- 
terns that vary in time, and phenomena caused by interfer- 
ence we can consider two fixed patterns that have different 
pattern speeds. We use a spiral argument for the i-ih spiral 
density wave (similar to equation [7| 



where we have dropped the constant and the angle 8i de- 
scribes an angular offset at time t = 0. With two logarithmic 
spiral perturbations the mass surface density as a function 
of time 



£(r,0,t) = Ai(r)cosf>i) + A 2 (r) cos(<^>) 



(10) 



The amplitude functions Ai(r) we assume are slowly varying 
with radius. 

Consider a location and time where there is a local den- 
sity peak so that E has an extremum with ^(ro, #o, to) = 0. 



By taking the derivative of equation ( 10 1 the extremum oc- 
curs where 

Ai(ro)misin(<^i(r ,6>o,io) + A 2 (r )m 2 sin(^ 2 (fo, Oo, to)) = 0. 

We consider small perturbations in time and angle about 
this moment and position; i = to + dt and 9 — 9o + d9 and 
require that this new position and time is also a local density 
peak and so is also an extremum so that ^ (ro, 9o + d0, to + 
dt) — 0. By expanding this partial derivative to first order in 
dt and d9 we can estimate d9/dt. To first order the density 
peak moves with angular rotation rate 

d9 mi Ai cos 0if2i + m 2 A 2 cos 02^2 

~dt ~ 



mi A\ cos 4>i + m 2 A 2 cos (j> 2 

TUiAiQi + m 2 A 2 Q, 2 



(11) 



mi Ai + m 2 A 2 

where in the second step we have assumed that the local 
density peak is near the density maximum of both waves 
where cj>i ~ <f> 2 ~ 0. 

The above expression describes the angular rotation 
rate of the density peak. We see that its angular rotation 
rate lies between the the pattern speeds of the two density 
waves, with value depending upon the amplitudes of each 
wave. Consider a situation with the first wave faster than 
the second one, f2i > Q 2 , and with peak amplitude at a 
smaller radius than the second one. In this case Ai and A 2 
are smoothly varying with radius and Ai is high at small 
radii and A 2 high at larger radii. At an inner radius where 



the first wave dominates, equation (111 implies that f2 p is 
near fii but at a radius where the second pattern dominates 
fl p would be near fi 2 . The peak would effectively vary in 
pattern speed with pattern speed decreasing as a function 
of increasing radius. Density peaks with an angular rotation 
rate that depend on radius have been identified in simu- 
lations, for example see Figure 5 by Grand et al. (201! 
and Figure 8 by Grand et al. ( 2012b I. The maximum an- 



gular rotation rate is that of the faster pattern, and this 
is consistent with the angular rotation rate of the feature 



identified by |Grand et al. (2012b I approaching that of the 
bar at small radius (see their Figure 8). The variation in f2 p 
with radius causes the winding angle of the density feature 



to vary in time or wind up (as shown in figure 7 by Grand 
|et al.|[2012b| ). We note that other types of transient models 
might also predict such behavior, for example swing amplifi- 
cation models predict that the winding angle varies in time 
as the amplitude varies. A swing amplified model would also 
effectively give a peak with angular rotation rate that varies 
with radius. 

The above expression for the angular rotation rate 



4>i{r, 9, t) = ai lnr — m, ( 



flit — 9i 



(9) 



(equation 111 does not depend upon at or the winding an- 
gles of the patterns. Consequently it can be used for a, = 
and so can describe a bar-like perturbation that has a fixed 
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angle as a function of radius. We can consider interference 
between a fast bar perturbation and a slow local spiral per- 
turbation. Outside a bar's corotation radius, the bar pat- 
tern speed would be faster than the local angular rotation 
rate. Spiral patterns are often observed in simulations to be 
strong inside their corotation radius (or at radii where the 
pattern speed is larger than the local angular rotation rate; 
e.g., Quillen et al. 20111. In a setting where a local spiral 



Here larger /3i , P2 correspond to more tightly wound struc- 
tures. 



pattern interferes with perturbations from a bar, we would 
expect that the effective angular rotation rate d9/dt of an 
interference peak between would lie between the bar's and 
the spiral's pattern and ranging at small radius with an an- 
gular rotation near that of the bar to near that of the spiral 
at larger radius. In this setting the effective angular rotation 
rate of the interference peak would necessarily pass a point 
where d8/dt is equal to the local angular rotation rate and 
so the interference peak would act as if it were going through 
its corotation resonance. This is likely the situation in our 
simulation and causing migration. Bar corotation regions are 
seen as location where migration is common even after the 
formation of the bar (in our Figure [I] but also see Figure 7 
by Minchev et al.|2012 \ . However slow spiral structure in the 
vicinity of a bar would cause interference peaks with nearly 
corotating angular rotation rates. This may in part account 
for the efficacy of bar corotation regions in causing stellar 
migration. 



3.2 Radial position of the maximum density 

Above we considered the azimuthal location of peaks. We 
now consider the radial location of the maximum density in 
a local peak and how it would vary in time. In our simulation 
we see local peaks usually appear at a small radius and move 
outwards, but sometimes appear at a large radius and move 
inwards. We can try to understand the speed and direction of 
motion using the same model as we used above. A localized 
burst of star formation could be caused from density wave in- 
terference. As long as the amplitudes are only slowly varying 
in time and radius then the local density primarily depends 
on the arguments. When the two patterns lie directly on top 
of one another then the density is a maximum. Consider a 
radius, angle and time where the two waves constructively 
add or 0i(ro, #0, to) = 4>2(ro,8o,to) = 0. Consider a nearby 
radius ro + dr where the two waves also constructively add 
but at time to + dt so that 0i()"o + dr, 80 + dO, to + dt) = 
4>2{ro + dr, 80 + d8, to + dt) = 0. We can solve for 80 + d8 
using equation (J9J) for both angles, 

80 + d8 = — i- ln(r + dr) + f2i (t + dt) 



mi 

0.1 

rn 2 



ln(r + dr) + n 2 (t + dt). 



(12) 



Expanding to first order and relating dr to dt, we find that 
the density maximum moves inwards or outwards in radius 
at a speed 



dr 



ro{^2 — fii) _ ro(n 2 — f2i) 



peak ai/mi— a 2 /m 2 ft — ,82 
where we have defined a parameter 



Hi 



= cot7i. 



(13) 



(14) 



We can see from equation ( 13 1 that the peak maximum 



moves slowly if the two pattern speeds are similar, and this 
is expected as they would remain in phase for longer. If 
we adopt Q.\ > Q.2 then the direction of motion depends 
on the sign of the denominator and this depends on the 
winding angles of each pattern (equation 181. If the faster 
pattern is more tightly wound than the other (/3i > P2), 
then the density maximum would move inwards, otherwise 
it would move outwards. In the first case we would expect a 
burst of star formation that progressively moves inwards in 
radius, and in the second case moving outwards in radius. 
Star bursts propagating inwards would occur if the center 
of the galaxy had more open spiral structure than the outer 
regions of the galaxy. 

The denominator of equation ( |13[ ) is only zero when the 
two waves have the same winding angle and so the two pat- 
terns are in phase at all radii simultaneously. If two patterns 
have nearly the same winding angles then the denominator 
is large. This implies that the peak is nearly in phase at all 
radii simultaneously. We would expect a burst of star for- 
mation across the density peak that is nearly coeval at all 
radii. We note that the above expression neglected variation 
of the amplitudes A\,A 2 with radius and time, but expect 
these would not strongly affect the drift rate, r, as long as 
the density variations are more strongly dependent on angle 
than radius. 

It may be useful to put equation (|13[) in units of pc/Myr 



'peak 



= 200 pc Myr 



-1 (fo-Qi) 
fi (/3i - P2) 



( ) 

V 200km s~V 



(15) 



There may be a progression of ages of star clusters and 
moving groups with stars closer to the Galactic center be- 
ing younger in the Sco Cen region (Eric Mamajek, private 
communication). Such transient peaks would be expected 
within an interference model if the outer galaxy is more 
tightly wound than the inner galaxy. 



4 MIGRATION BY LOCAL DENSITY PEAKS 

Above we have illustrated how interfering spiral density pat- 
terns can cause a local density enhancement to appear for a 
short period of time. We now consider whether such a peak 
can induce stars to migrate. We first review a Hamiltonian 
model for migration caused by a corotation resonance based 



on the seminal work by Sellwood & Binney ( 2002 1 



4.1 Corotation resonance model for migration 

For an axisymmetric system we can describe orbits in terms 
of a Hamiltonian that is a function of the angular momen- 
tum, L, and the momentum associated with epicyclic mo- 
tion, J. These two momenta are conjugate to the an angle 8 
that is approximately the azimuthal angle in the plane, and 
the epicyclic angle, ip; 

H (L, J; 8, <p) = g (L, J) ~ g (L) + 9l (L)J + g 2 {L) J ^. (16) 

As the system is axisymmetric, the Hamiltonian does not 
depend on the two angles, so Ho, J and L are all constants 
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of motion. Above and on the right we have expanded the 
unperturbed Hamiltonian to second order in J. The angular 
rotation rate of a particle in a circular orbit in the unper- 
turbed potential, Sl(L), is 



fl(L) = 



dgo(L) 
dL ' 



(17) 



The epicyclic frequency «(L) = gi(L). For a flat rotation 

curve g (L) = v 2 lnL, Q(L) = and k(L) = V2Q(L). 

We can consider a logarithmic perturbation to the grav- 
itational potential in the form 



V(r,0,t) = e/O8Inr-(0-n.t)) 



(18) 



where r is the radius, the perturbation has amplitude e(t), 
pattern speed Sl„, and winding angle 7, with cot 7 = p. We 
consider a system with 9 increasing in the direction of ro- 
tation. In this case a trailing logarithmic spiral has j3 < 0. 
The above perturbation depends on radius rather than our 
action angle variables. If we take the low eccentricity limit 
and average over the dependence of r on the epicyclic angle 
then r ~ L/v c where v c is the circular velocity with angu- 
lar momentum L. In the vicinity of a corotation resonance 
(and distant from any Lindblad resonances and in the low 
eccentricity limit) the Hamiltonian depends on angular mo- 
mentum L and azimuthal angle 9 alone; 

H(L, 9) = g (L) + e(t)/(/31nL - (9 - O.t)). (19) 

where we have replaced r with L and removed an arbitrary 



constants from 9 and t. Sellwood & Binney (20021 have il- 



lustrated that in the vicinity of a corotation resonance, the 
epicyclic amplitude is not affected. 
It is useful to define an angle 



<t>= p\nL-{e-Q, s t). 
with time derivative 

4> = Py -(6 -Sis). 



(20) 



(21) 



The potential perturbation function, f(<fi), may be periodic. 
For example the function f(<j>) = cos(m0) gives an m armed 
spiral structure. This type of perturbation is equivalent to 



the model discussed by Sellwood & Binney ( 2002 ) and shown 
in their Figure 6. Previously we discussed density perturba- 
tions using a similar form for the perturbations. Here we 
focus on perturbations to the gravitational potential. 

Hamilton's equations (using the Hamiltonian of equa- 



tion 19 1 are 



L = 



dH 
~~d9 
dH 



Inserting 9 and L into the expression for tj> gives 



n(L). 



(22) 
(23) 

(24) 



The above equation implies that when fl(L) ~ fl s (in the 
vicinity of corotation), the argument <f> is nearly constant 
and the angular momentum increases or decreases depend- 
ing upon the sign on the right hand side of equation ( |22[ ). As 
ef(4>) is the perturbation to the gravitational potential we 
can also write equation (22 I (one of Hamilton's equations) 
in the following way 



t = ~ = -(rx V#) z (25) 

where $ is the entire gravitational potential. The angular 
momentum drift, L, arises from the torque from the non- 
axisymmetric portion of the potential perturbation. Lagging 
the potential minimum we expect e/'(0) > corresponding 
to L > 0. Particles lagging a potential minimum would mi- 
grate outwards. 

A maximum migration rate is 

L max ~ max \ef'(<j>)\ (26) 

leading to a maximum radial migration rate of order 

. L max 

1 max • I ^ 1 I 

Vc 

Note that the maximum torque on a particle can be directly 
estimated from the maximum tangential force at a given 
radius. 

In the rotating frame with 9' = 6 — Q B t, the Hamiltonian 
becomes 



K(L,9') = g (L) - LSI, + e(t)f(/3\nL- 9'). 



(28) 



As K is independent of time, it is conserved and known 
as the Jacobi integral. Expanding this Hamiltonian about 
the angular momentum corresponding to corotation or L — 
L s + J, with Q(L B ) — Q, s , 

T 2 

K(J,9') = const - fi'(L s )y +e/(/31nL s -9 ) + .... (29) 

with fi'(L) = —v 2 /L 2 for a fiat rotation curve. This gives an 
equivalent description of the migration model in the coro- 
tating frame. Up to this point we have not assumed any 
particular form for the potential perturbation, only a pat- 
tern speed. Nevertheless the maximum migration rate can 
be estimated from the maximum strength of the derivative 
of the potential perturbation (equations 26 271. 



4.2 Corotation resonant width and eccentricity 
limit 

In the WKB approximation a logarithmic spiral density per- 
turbation comprised of a single Fourier component 

E(r, 9, t) = -SVnEo cos(a In r - m(9 - O s t)) (30) 

has gravitational potential 

V(r,9,t) w 2itGS ™ Eo cos(a\nr-m(6-n s t)), 
\a\r 

(computed with wavevector k ~ a/r as by Binney fc 
|Tremaine| 1987). Here 5* m So is the strength of the m-th 
Fourier coefficient of the surface density. Eo is the az- 
imuthally averaged density and S m is an amplitude in 
units of Eo. Previously we have used a winding parame- 
ter (5 = a/m so that the argument can be written a In r — 
m(9 — fl„t) = m[p\nr — (9 — fl s t)] — mcj). In this case the 
function / of Hamiltonian equation (1191) is a cosine, and the 



(31) 



perturbation strength for the potential (equation 18 1 
27rGS m E r 



(32) 



The Hamiltonian (equation 29 1 is similar to that of a 
pendulum. Shifting the angle 9' by a constant we can rewrite 
the Hamiltonian as 
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K(J,6') = 



L'i 2 



i{m0') 



(33) 



where we have neglected the constant term in the Hamil- 
tonian, assumed a flat rotation curve and flipped the sign 
from equation ( |29[ ). The resonance width can be estimated 
by computing the energy of the separatrix (K = e at J = 0, 
m9' = 7r) and then solving for J at 9' = 0. The width of the 
resonance (peak to peak in J) is 



(34) 



Here A J describes the width of the corotation resonance and 
determines the maximum radial distance that a particle can 
migrate in the resonance; 




(35) 



where r a = L 3 /v c . Using equations (321 and (35 1 



Ar 



0.9kpc ( 

(£)"'( 



015 



E 



8kpc I V50M Q pc 2 



(36) 



-200km s" 1 , 

Here we have used a surface density similar to that of the 



Milky Way disk at the Sun's galactocentric radius (Holm- 
berg & Flynn 2004), a parameter a = 17 for m = 4 and a 
winding angle of 13° as preferred by some models for the 
Milky Way ([Vallee |2008[), and the 15% density contrast es- 
timated by |Drimmel fc Spergel| ( |2001[ ). This suggests that 
unless the Galaxy experienced strong and open spiral arms, 
migration over large distances is unlikely. 

Using the cosine perturbation we can also estimate the 
maximum migration rate using equations ( 26|27 1 and the 
derivative with respect to 9' of equation (31 1, 



2-KGS m T, r s 

P 

2km 



i / SV^\ / J^_\ / So 
V0.15/ ^8kpc/ I^OMqpc 2 

\4/ V 200km s _1 



(37) 



(here a = m/3). Only if there are strong and open spiral arms 
is migration rapid. We see migration at a rate of about 20 
km s _1 (see Figure [5]) exceeding the value above. Grand et 
al. (2012b I show a migration rate (using their Figure 10) 
of about lkpc in 20 Myr equivalent to 50 km s _1 . The mi- 
gration seen by Grand et al. ( 2012a|b| and us is both more 
extreme and quicker than estimated above using parameters 
estimated for local Galactic spiral structure. We will discuss 
this issue in more detail below. 

For a particle to migrate the distance given in equa- 



tion (35 1, it must be captured as the pattern grows near 
<f> ~ — 7r/2 and at low L and then released at <j> ~ 7r/2 at 
higher L, as described by Sellwood & Binney (2002). The 



dependence of the resonant angle on /3 In L does not affect 
the resonance width estimate (one can follow the same pro- 
cedure to estimate the width of the separatrix). However the 
shape and normalization of the function f((j>) does affect the 
resonance width. 



How low an epicycle (eccentricity) is required for cap- 
ture into a corotation resonance? The eccentricity must be 
low enough that the angle, 4>, remains within 7r/2 of the 
potential gradient extremum, otherwise the torque on the 
particle would oscillate instead of remaining positive or neg- 
ative, as exhibited by the migrating and non-migrating par- 
ticles shown in Figure |4j We consider cj> as a function of r 
rather than angular momentum or <f> — /31nr — 9'. An or- 
bit with eccentricity e has an apocenter of r ~ ro(l + e). 
Inserting the apocenter radius into the expression to <f> we 
find (j> ~ P lnro + Pe, consequently epicyclic variations cause 
a change in (j> of order 0e. To keep <j> from varying by more 
than 7r/2 as the particle movies radially in the epicycle we 
require that 



7T 

Tp 



ixm 
2ctT' 



(38) 



The corotation resonance is likely to be ineffective for par- 
ticle eccentricities above this value. Tightly wound patterns 
can only cause migration in extremely low eccentricity stars. 

Previous studies have found that corotation resonances 
are only effective at causing migration among low eccen- 



tricity particles (see Figure 12 by Minchev et al.|[2012 and 
associated discussion). Using the ai = 6 parameter for the 
m — 3 pattern from Table [T] listing our model patterns, 
equation (38 1 gives a maximum eccentricity of 0.8 for par- 



ticles to be captured into a corotation resonance. However 
in section [2~Tj we estimated that only particles with eccen- 
tricity e < 0.2 migrated, hence equation (38 1 overestimates 



the maximum particle eccentricity. We will discuss this issue 
again when we consider the Gaussian bar potential pertur- 
bation model below. 



4.3 Migration by a thin Gaussian bar 

When the potential perturbation is well approximated with 



a single Fourier component (equation 31 30 1, then the max- 
imum migration rate occurs at an angle <j> — tt/2. However 
the density as a function of angle is poorly approximated 
by single Fourier component (see Figure [3| and we find that 
migration is taking place when the particle is 30° (not 90°) 
of the density peak. Similarly the extreme migrators shown 
by Grand et al. (2012b I were within 10 - 20° of a density 



peak (see their Figure 14). We could expand the gravita- 
tional potential perturbation function in terms of Fourier 
coefficients with f((j>) = A m cos mc/>. However the maxi- 
mum tangential force (and the angle at which the maximum 
is located) is then a function of the Fourier coefficients as 
the derivative of each component reaches a maximum at a 
different angle. 

To take into account the narrowness of the spiral fea- 
tures, we describe a local density peak as a bar-like linear 
feature with a surface density profile (rather than in terms 
of Fourier components) 



-y 2 /(.2<? 2 ) 



(39) 



where E p is the peak density subtracted by the mean density, 
and the density bar has full-width-half-max of 2.35 a. This 
description has the advantage that it is a function of only 
three parameters, the peak density, E P) (above a background 
level), a width, described with a, and an orientation angle. 
We assume that the density feature does not strongly depend 
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on the direction perpendicular to y. We orient the density 
feature with the winding angle 7, with 7 = corresponding 
to the feature azimuthally oriented and 7 = 90° with the 
feature oriented radially. Here y increases in the direction 
perpendicular to the feature. For a tightly wound structure 
P^ 1 = tan7 ~ 7. 

The Fourier transform of the above density profile is also 
Gaussian with E(fc) = E p <7 exp(— k 2 a 2 /2). In the thin disk 
approximation, and using Poisson's equation, each Fourier 
coefficient of the gravitational potential tx e !fcH_ ' fcz '. Using 
a pillbox about the plane we can show that for each Fourier 
coefficient the gravitational potentiaQ 



W \k\ 



exp 



1 2 2 
-K a 



+ iky 



Taking the y derivative and integrating over k 



9$ , — 

= 2v27rGE„ 

dy 



-k 2 /2 



sm ( — j dl 



= 4^/%GTjp DawsonF 



V2a 



where the Dawson-F integral is defined as 



F(x) 



exp(— x 



expy dy. 



(40) 

(41) 
(42) 

(43) 



The Dawson integral peaks approximately at F(l) ~ 0.5, 
consequently the maximum drift rate should occur for a 
particle located at a distance x ~ %/2cr from the peak or 
approximately at at the half-width-half max position. 

We can now consider a rotating logarithmic Gaussian 
bar with surface density 



E(r,9,t) = E p exp 



Cfllnr- (6-nt)) 2 
2s 2 



(44) 



with angular width s. We can relate s — a jr. For a Gaussian 
bar tilted at an angle 7 we can approximate locally 

9$ <9<3> . 

The angular derivative of the gravitational potential 
9$ 



(45) 



-' 4 v / 7rGE p sin 7 DawsonF . _ 
of) V v2s 



(46) 



with 7 the winding angle and <f> defined as in equation ( 20 1 



Using the derivative d&/d6, we can estimate the max- 
imum torque T max and from this the maximum drift rate 

r ~ T m ax I Vc OV 

~ GE p 2 v / ^~ 1 sin7 (47) 



or 



12km s~ 



(skpc) ( 



50M Q pc- 2 

Vc 

200km s" 



/ sin 7 \ 
V 24° J 



(48) 



Here we have used the typical 24° winding angle estimated 
for the spiral features in the simulation ( Quillen et al.|20~TT I 



1 The integral of <3?(fc)dfc does not converge, however the integral 
should be cut off at small k due to the finite thickness of the disk. 



and a surface density similar to that of the Milky Way disk at 



the Sun's galactocentric radius (Holmberg & Flynn 20041 



The above estimate suggests that high density peaks are 
most effective at causing migration. The migration rate is 
only fast enough to cause significant migration within an 
orbital period if the surface density is of order 50M© pc~ 2 . 
As mentioned previously, this presents a potential problem 
for migration models for the Milky Way as peak spiral am- 
plitudes have been estimated at levels much below this (e.g., 



Drimmel fc Spergel|200l| ). 

The maximum drift rate occurs for a particle approxi- 
mately at the half-width-half-max location and so is sensi- 
tive to the peak width. The particle must remain within a 
half-width of this location, hence the upper limit on eccen- 
tricity for migration particles must be modified from that 
estimated using a cosine potential perturbation (equation 
381. Using an apocenter radius r = ro(l + e) and keeping 
<j> within angular width s of the maximum migration rate 
position we estimate that for a particle to migrate 



(49) 



For fi ~ 2 in our model (taking values from Table [T]) and 
width s ~ 30° « 0.5 radians from Figure [3] we estimate e < 
0.25 and this is approximately consistent with our previous 
estimate (based on Figure [2^,) that the maximum migrator 
eccentricity e max ~ 0.2. 



Our estimate for the migration rate (equation 48 1 de- 
pends on the strength of the potential term and this is set by 
the peak surface density of the feature. Using the Gaussian 
bar model we compare estimated migration rates to those 
measured in the simulation. In Figure [5] we show both the 
migration rate and density of the nearest spiral peak for 
our selected particle (with trajectory also illustrated in Fig- 
ure[2^, [3] and[4|. As expected, the migration rate is tightly 
correlated with the peak density of the nearest spiral fea- 
ture. During the period of migration from 500 Myr to 700 
Myr, the density E p of the nearest density peak has val- 



ues in the range of lOO-25OM pc" 2 . Using a mean value of 
E p w 18OA/0 pc -2 and equation (48 1 we estimate a migra- 
tion rate of 40 km s _1 . This is somewhat higher than the 
mean migration of about 20-30 km s _1 , but we should re- 



member that equation ( 48 1 estimates a maximum value and 



so we expect the actual value to be somewhat lower. 

We have found that the Gaussian bar model success- 
fully estimates the migration rate, the limiting eccentricity 
value for migration and the angular position of a particle 
that is rapidly migrating. In these ways the Gaussian bar 
model is superior to a model that approximates the potential 
perturbation with a a single Fourier component. However, 
even though the derivative of the potential perturbation is 
well defined for the Gaussian bar, the potential itself is not, 



making it impossible to modify equation ( 35 1 to estimate 



the extent that a particle can migrate. Up to this point we 
have not considered how growth and variations in winding 
angle effect the migration model. However if migration is 
rapid, then a local approximation for the Gaussian bar is 
sufficient to estimate the migration. To estimate the extent 
that a particle can migrate we would have to take into ac- 
count the time dependent form of the potential perturbation 
and the particle's trajectory (in <j>) with time. 
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Figure 5. Peak density and migration rate. For the selected 
particle shown in Figure [4^, the nearest peak mass surface den- 
sity density, E p , subtracted by the mean density, is shown as a 
red line and has y-axis on the left. Also shown as a green line and 
with along with j/-axis on the right, is the rate of radial migra- 
tion, r. Both S p and r are shown as a function of time. During 
the period of migration from 500 to 700 Myr, the peak density 
and migration rate are correlated, as we would expect from a 
corotation resonance mediated migration model. The rate of mi- 
gration and peak density are approximately consistent with the 
maximum estimated using the Gaussian bar model (equation|48|. 



5 DISCUSSION AND SUMMARY 



As have other studies ( Grand et al.|2012a|b| |Minchev et al 
20121, we find stars that migrate in radius in our N-body 



simulation. However spectrograms of the simulation do not 
show transient patterns that are strong near their corota- 
tion radius, presenting a difficulty for the corotation reso- 
nance model for mig ration ( jSellwood fc Binney||2002} . We 
find that particles migrate outwards in the simulation when 
they are at low eccentricity and lag a short lived local den- 
sity peak, confirming the results of previous studies ( |Grand| 
et al.|2012a|b l. Density peaks may not survive more than a 
rotation period, suggesting that they are not due to individ- 
ual growing and decaying spiral density waves. We consider 
the possibility that they are due to interference between den- 
sity waves. This scenario does not require that spiral density 
waves last with fixed amplitudes, pattern speeds or winding 
angles for dozens of rotation periods and so is not inconsis- 
tent with recent studies that suggest that spiral structure is 
transient (e.g., Sellwood 20111. 

We have explored a simple model for interfering spiral 
density waves and find that they can produce short lived 
local density peaks and these, if they are sufficiently strong, 
can account for the migration seen in our simulation. The 
mechanism for migration due to a corotation resonance pro- 



posed by Sellwood & Binney (2002) can be modified to apply 



for a short-lived local density enhancement. Using a Gaus- 
sian bar model for the potential perturbation, estimates of 
the migration rate, angular offset between particle and spi- 
ral feature, and maximum eccentricity for migrators roughly 
agrees with the values measured in the simulation. We show 
that interference peaks can appear to rotate faster at smaller 
radii than outer radii (and so appear to wind up) as seen in 



our and other simulations ( Grand et al.|[2012a |. The angu- 
lar rotation rate of an interference peak should lie between 
the pattern speeds of the two patterns that cause it. Nearly 
corotating interference peaks would be particularly likely to 
appear near the end of a bar due to interference between the 
fast bar and more slowly moving spiral patterns outside the 
bar. This may account for ongoing radial migration that is 
seen in simulations after bar formation and in the vicinity 
of the bar. 

We can contrast a migration mechanism mediated by 
individual spiral density waves (as proposed by [Sellwood fc| 
Binney||2002"t and one mediated by short lived interference 
peaks. An interference peak can have a higher surface den- 
sity than each individual pattern. As the distance migrated 
and speed of migration depends on the surface density, in- 
terference peaks can enhance migration. Many interference 
peaks can be produced over the lifetime of a few spiral den- 
sity waves. This suggests that migration events could occur 
faster and more frequently than when considering a model 
where migration events only occur when a single pattern 
grows and decays. This may ameliorate the timescale prob- 
lem as only a few dozen spiral density waves could arise and 
decay during the age of the Galactic disk, however a few 
times this number of interference peaks may arise during the 
same time if there are multiple spiral density waves present. 
With a single pattern, migration can only occur at the coro- 
tation resonance of the pattern. Interference peaks can have 
a variety of effective angular rotation rates, depending upon 
the amplitudes and patterns of the waves causing them. Con- 
sequently migration induced by interference peaks may be 
more pervasive than migration induced by individual spiral 
density waves. 



The migration rate and the distance a star can migrate 
due to corotation resonance are dependent on the peak den- 
sity in a spiral wave. As the stellar surface density decreases 
with increasing radius, migration due to corotation reso- 
nances may be less effective in the outer parts of the Galaxy. 
Furthermore as migration likely occurred near the Sun's ra- 
dius it is likely that spiral structure is either stronger than 
previously estimated (e.g., by Drimmel & Spergel 20011, 



or was stronger in the past than currently observed in the 
Galaxy. Better measurements of the spiral structure in the 
Galaxy are needed to make quantitative estimates of the 
extent of on-going radial migration. 

Interference between spiral density waves can mimic 
transient-like behavior, including inducing migration as dis- 
cussed here. Short lived density peaks can also cause local- 
ized bursts of star formation that can move across a galaxy. 
This picture of star formation differs from that expected 
with a single strong density wave that would produce a con- 
stant amplitude wave of star formation. We estimate that 
the rate that a burst of star formation moves across a galaxy 
depends on the winding angle of the two patterns. A burst 
would move inwards if the outer galaxy contains more tightly 
wound spiral structure than the inner galaxy. Future studies 
could probe for evidence of localized bursts of star formation 
in context with spiral structure models to better constrain 
the nature of spiral structure in galaxies. Future studies can 
also aim to quanitify the statistics of migrating stars using 
interference mediated corotation migration models. 
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